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ARSTRACT 


A finite element method Tor structural opuimigevis oes 


Mepicitavec bieam ol a Nnomogeneous, 1SOtCrOopiec maver iia. 


tr) 


developed. The beam has a rectangular cross-section cz 


menstant fixed height, a fixed length, and a fixed volume. 


tf) 
9 


Seructural optimum is defined as that shape which allcw 
maximum load within the elastic range. 

PecOmOuGerm Drog@ran as cdevelooned to solve the resus ine 
Bpyevem Of equations and various example problems are s:ived. 
Comparison is made with exact optimum beam designs wheres 
possible. 

The finite element model is able to solve problems with 
aay boundary conditions and types of loading that are zon- 


Sistent with the number of elements selected. 
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Le INR UC IIa: 


Optimization has a history as old as the universe itself. 
Whe planets have taken their optimum position in the solar 
System, the solar system in the galaxy, etc. The whole 
Mmesvory Of evolution is one of gradual optimization. Man 
has used optimization in one form or another since the 
earliest of times. Individual man was weak with respect to 
Mis hostile environment. He learned to become a social 
animal so that as a group he was strong - optimization. 

In wars between groups spears defeated rocks, arrows defeated 
Ties 2G, SO Onavo the preseny day. Alle were Casenuaa a 
SovamaZation. 

The early Egyptians formalized some optimization tech- 
niques in their development of geometry and trigonometry. 
We development of the Calculus enabled us Go find optimum 
momel Oneawedmmes. ThemGalculus of@Varietaiens Magdemuecs! bie 
mite, ODLIMmiZatvion of functionals, the optimum funetion from 
a function space. 

iMipsmicve LUica @OouIl Zol1@n, Mewever {Ol Tenvimes Une 
falculus of Variations produces nonlinear differenvial 
equations that cannot be solved in closed form and the 
mumer1 cal methods of solution are difficult. if these 
differential equations could be transformed into algebraic 


equations, a solution might be more readily obtained. 








The finite element method has been arvlied to the 
theories oi elasticity {lj}, [21] and mechanics io 
[3], as well as to other fields, to resolve such problems. 
mc has been used ian the optimization of truss neg oer ogee 
as yet not in actual beam optimization. Since the method 
1s essentially variational in nature, it should have an 
peeericavion to structural optimization theory. This @suc— 
Bests the desirability for research in this field. 

pimce Ciiceisewa first attemply at ae tormularrcnmemerr jane 
ype, we will restrict ourselves to a specific type of 
mecowp lem; the structural optimization of a prismatic beam 
@f One material which is homogeneous and isotropic. This 
beam has a rectangular cross-section of uniform fixed height, 
and has a given volume and length. The structural optimum 
in this case is the shape required to maximize the load 
within the elastic range. 

A variational formulation is first presented ta estaplisc, 
Mrewoptimigation technique for isoperimetrre problems. The 
finite element method is applied to the problem using a 
parallel approach. A model is developed and then checked 


for validity by solving various applicable problems. 


18, 








Ti. VARIATIONAL FPORMULATION® OF BEAM OPDRUMTZ ier Orie iis 


foe DER INIREON OF THE PROBLEM 

Salinas 44) Shows that the probllemeor seruce une ieee ce 
MiZavion Of a beam, under a fixed volume Constrainc, may be 
solved by using calculus of variations techniques. He shews 
mea. OPeravion wWwiuh the Cechniques for 1soperimetric prep— 
lems [5] on the augmented potential energy funce ones 
T* = T-AS, produces equations which describe the struc- 
meralivy Optimum beam. T is the potential energy of the bean 
and S is the isoperimetric constraint. The quantity A is 
a constant Lagrange multiplier which is actually a measure 
of the strength of the beam. 

We shall consider the problem of Structural opuimZavien 
ef a simply supporcted beam of length, L, and volume, ie 
ie applied load is uniform, - (Porce/inie. lengcn sire 
cross-section of the beam is such that the moment of 
inertia may be given as I = CA" where C and n are constants. 
Gravity forces and shear effects are assumed negligible and 
mee usual beam assumptions 6) are used. 

i) small deflection theory 
11) Winlasmea ly ouress state 

Figure 1 graphically depicts the beam and defines the x 

and y axes. We seek the shape of the beam to give a maximum 


lead amcensi1ty subject to the constant volume censeraing. 


11 
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| | pent ea SES ! 4 SSE ETE pee ee r 


ys a ae, cee ea ae 
Lae | ! | | . 
2 ee ee ee eee 
ae. pm ae OTT <3 x 





A(x) (not necessarily rectanrular 


Papere ll. iosdmpbly Supported Bean 


B. THE VARIATIONAL FORMULATION 
Wemiarst Form the potential energy TPuncevonal iacem 
strength of materials considerations. The displacement of 
the beam is denoted by v and the primes indicate differen- 
PiewemOrl With respect to x. 
ie , PE co? = PY | dx 


0 


L 
iene d suetnimne scons traim@may peawri oben asae A(x)dx=V. 
0 


thus the constraint augmentation becomes: 


(6 
oe er A sy ed x 
@) 


Noting that I(x) may be denoted as CiG =. the augmented 


movential energy functional is: 


L n 
We fee | BEAG) Gate ee ay Aa) | ‘elb< 
0 c : | 


where v and v" are v(x) and v"(x) the displacement and 


curvature of the beam respectively. 


Ae 








In order to extremize the above functional, we first 
wake the variation. [7] of the functional wilh récnecrare 
displacement, v, and set this equal to zero. The result 
is the differential equation of equilibrium. This variation 


yields: 
n ! 
Cole eee) P= 0 
which may be written: 
1! 
(ATy") = P_/EC 
itinerrating this equation twice with respect to xX gives: 


2 
X 
" = —<— 
Vv P i EX 5 +- Ky x +- K,, 


ae 


MjeSweonstcants Of Integration. Ky and Ky > are obtained 


meom Lhe boundary conditions. For a simply supported beam, 


the end moments are zero, or: 


(a) A’(O)v"(0) = 0 
(b) A’(L)v"(L) = 0 
Boundary conditions (a) and (b) yield K, = 0, K. = -L/2 


2 


respectively. 
Thus we obtain the differential equation of equilibrium 


mole une Simply Supported beam. 


|e 
any" = <8 Cage Ghee 


Taking the variation of the funetional wren respecu me 
deco x ye and setuine this equal to Zero, 2ives tiewcia 


ferential equation which we refer to as the optimality 


iS 








condition, For, our exarole the ontiralitysconc trae ee 


NEG cayD-lo yn 2 sais 


5 = (0) 
So line ene Sy ie lde : 
- (+) 
ft 2 2 
ae Or A (TI-2) 
nec 
n-l 
espe 2. a 
A We = Vel Cit =2e5) 


Equation (II-2a) shows that, for an optimum beam, the 
eurvacvure 1S constant only for n = 1. Otherwise the producer 
®t a power of the cross-sectional area times the curvature 
imemeonstant . 

Multiplying Boeiwsot des of equation Cit= 2) by A’ and sub- 


stituting into equation (II-1), the development continues as: 


oad 


nt] 
Nowy 2 e = 
Avv" = A nEC CL = 
From equation (II-1): 


ideal 


|e 
2 On O 2 
A” Jagex = deg (6 Lx) 


Hence , 
2 
ale: 
alae A O a 2 
eee ee 
; 2 ae oe OE aa : ‘ ; 
Since (ee =. = CL xx) Pit 1s Convent tent tame ere 


development to use the latter expression. Substituting this 


Mite chne previous equation yeilds: 


yh 











Taking the (n+1)°A root of both sides 


nPé 2 Grimes 
A = cael CE x=x 5) (ia 
Introducing the constraint equation and solving for iA we 
obtain: 
5 i 
L a nP 5 a Meee ik 
: A(x)dx = 4 BEC Clee ) a= ie 
Bue since A is a constant we obtain, 
i 
—— = “ 
B@ar ak Ee. Se 
d uk: 1 \ pe 2s (II-5) 
a OFC We Ol x = ) ox 
010 
Rearranging (II-5) and solving for er we obtain the 
equation defining the maximum load intensity 
ntl 
r 
=a ary II-6) 
ae 7 7 iz igaede Gis 
ip (Lx-x* an ax é 


0 


Siiaebitcucine (8t=-5) intow(fi—4) we cbtainm for tae opel 
mum area: 


(le Tea) 
f (ox-xe "tla | 


i> 








It should be noted here that the interral is nevada: ee 
in closed form for all n and therefore must be numerically 
approximated for some problems. 

Summarizing the equations that define the solution to a 


given problem we have: 


2 
Vo (Lx=x® ) 9"? 
a ai al MESS) 
i: a 
df Cine” nt lay 
@ 
ae eal 
ECR (y ) 2 
oh : 2 ntl GE) ) 
Ji Cuca ntl a, e 
0 


The value of A is determined as follows. From equation 


(TI-2a): 


J nEcan-t 


PeOUecmECietmnol MatverlLals= eneory : 


M eli ) 
tee. max 


If we define the yield stress of the beam material, Sy as 


the maximum allowable stress, we obtain: 


where See Ve, Diner vertical Gistance to the tieunesly ibe 


from we meutral axis. Thus: 


116 











Ei ( Dg) eee 
we max max © pyn - 
VY T max max 
Sibstitutane for’! 
2 in ae a 
J ECA” 
max 
© oom | ae 
n n-l 
y CA max 
ns* cannl 
a (Timon 
OF Re 
max 


n-l 
It may be shown by example that the quantity a. es 





mmaependent of x and thus that A is indeed a constant for a 
Given type of cross-section. 

The first example is a Soren height rectangular cross- 
Beevion Of variable width. The height and width are h, and 


ibex) respectively. 


= O 
16 2 pin” aca “Bes 


Lee xh 


o 
Te - A(x) 


AL 








Thus for a rectangular cross-section as e2bcevenseouacued 


oli—-10) becomes; 


NO 


S 
AN = ; (II-11) 


is 


ea) 


ine second example is°a- Circular cress—sectuonmes 


Pecgius r(x). 
ee 2,2 1 
I(x) = PEG 2 (ar(xy?y° + (4A) 


Rese) = Tees 





Haus for a circular cross-section, equation (11-10 


becomes: 


. 


(Tse 


J aay 


E, 


The third example is a rectangular cross-section of con- 


stant width, b, and variable height, h(x). 


3 
Dig x 3 3 al 
iG) Se eh 
ie ate 
Mex. )ea=- Dome) 
Lee ts A(x)? 
12b 


ie 














C=—., n=3, R=h/? 
12 
oabe> Mace. ee en?) + {4)e2 
2 2 2 ia 2 
R R 12b h 


Thus for a rectangular cross-section of constant width 


ama variable height, equation (11-10) becomes: 


je eee (iii) 


Note that equations (II-11), (II-12), (1I-13) depend 
ary On cross section type and not on the manner of loadings 
(independent of the term r P(x)v(x)dx). They are also 
independent of the boundary conditions since they were 
wemeved [rom equations in which the boundary conditions had 
mee been envoked. This observation will be Used in Section 
my. 

Equations (II-8), (11-9), and (ITI-10), when applied to 
meee citic problem, yield the structurally optimum beam 
shape for maximum load intensity for beams which conform to 
the assumptions made at the beginning of this section. These 
beams must be simply supported and under a uniform load 
eeeecLribution. 

Appendix A gives three examples of structural optimiza- 


iG tena problems Sollee teem ey Ons ul — omc alee 


ey, 








III. A FINITE FLEMENT MODEL FOR BEAM OPTIMIZATION 


eo SPECIE CATIONS AND ASSUMPTIONS 

This development parallels that of the continuous vari- 
euronal formulation of the preceding seevion. We consider 
prismatic beams of one material which is homogeneous and 
meovropic. The modulus of elasticity and yield stress ef 
the material are E and Sy respectively. The cross-sectional 
area of the beam is rectangular with constant height, h, and 
variable width, b(x). The length and volume of the beam are 
denoted by L and pe respectively. 

The plane of loading is in the centroidal plane parallel 
Memene y axis. Gravity forces- and shear effects are assumed 
to be negligible. The usual beam assumptions [6] are taken: 

i) small deflection theory 


i.) mbes i irecer so bace 


feo ATEMENT? OF THE PROPLEM 

MomigsecllOnmll. We Seek tO Optimize the Shape ‘ol the 
beam (cross-sectional area as a function of position along 
the beam axis) to give a maximum load intensity subject to 
ee constraint that the beam volume is fixed. The finite 
element method results in an approximate solution to this 


problem. 


20 








GC. DEFINIVIO) OF Dew 1 oe ee 

The beam of Figure 1, without suppores, 15 discrerize: 
into finite elements of equal length as shown in Figure 2. 
The elements are identified by numbering them sequentially 
meom the orifin. The x (centroidal) axis and@yset- ew mocu 
tive downward) define the plane of loading and bending. 

We consider a typical element, say the pela as shown in 
Figure 3. The displacements and slopes at the left and 
ree hy ends of the element are denoted by Vi 28) > Vo > 8 
respectively. The cross-sectional area is A. ine woe aul 


coordinate system x. and y, are defined for each element. 


De APPLICATION OF THE FINITE ELEMENT METHOD 

In accordance with the finite element method [1], shape 
functions are assumed for the displacement of each element. 
memce a Cubic polynomial function for Vas the displacement 
over the ae element, yields continuity of displacements 
and slopes at each element boundary, it is taken as the 


assumed displacement field. Thus the displacement over the 
th 


at element may be written: 
v(x,) = <N(x,)> {v}, 
Loge eal 


where SNE is the cubic shape function vector as developed 
in Appendix B-1 and tv}, is the column vector (v,, 85> Vora 

The elemental cross-sectional area, Aes is assumed to be 
constant over the element to obtain a tractable mathematical 
Gece l. weEhis DOIN’ 1s further diseussed laver"in tine 


development. 


ae 











Figure 2. A Finite Element Beam. 


a 


Figure 3. The i Finite Flement. 


a 








Since the augmented potential enere sic ered cle 
T*® = U-W-AS, must be extremized, finite element formulations 
must be obtained for the: 
i) strain energy of bending (U) 
ii) potential energy of external force (W) 


iit) Wa sOperimetri C Ceonstralnve Gms) 


ieee cein Bnerey sof Bending 


The strain energy of bending of a beam may be given 
as; 


L 
Geof Se Gay oe: 
0 





2 


The finite element assumption that each element is 
under uniaxial stress enables the strain energy in the see 


element, Us; to be expressed as: 


A EI. 5 
uz. = f —— (vi')° dx. 
a 2 il al 

0 
where Ve denotes the second derivative of v(x, ) with respect 
ior X.. 

1 

Por a rectangular section of constant herehu bo and 


maeach , Ds > the moment of inertia of any section of the 


element may be written as: 


I. = CA, C = h°/12 
i i. 
The strain energy in the eee element is obtained as 
follows: 
= > 
v(x, ) <N(x, ) {v}, 


as 








! = aT ff * \ 
V (x, ) <M Once 


tf fi _ a Pt tt Bb 
(v ) =e ig <}y i 

A 
a no A eho 
ze 0 C 1 
1x4 


pineCeakWen produce 
shape function vector, <Ny 
Xs in the expression, the 
outside the integral, and 


inus . 


Jt 
A. {v}. / 


The integral is evaluated 


a 4Ux4 matrix. Therefore: 


_ EC 
Bi a Be 


Here it may be noted that, 


Meevor associated with it 


{fy} 
i 


meliues for the vector {v}i 


could not be removed from the integrand. 


ioe 
ab 
<N a tv}. 
tf A 1! 
<N SN Sy) oe 
lg a al al 
saa 1x4 Uxl 


of the second derivative of the 


2 N">,is the only function of 


remaining terms may be brought 


the integral may be evaluated. 


<Not 


<a aye ee 
al i aL i 


in Appendix RB-2 and called [k*], 


{v}; [k*] {v}, 
ibe ae Se 


Nx] 


if the A. has a shape function 
other than @ constanv, thesvec vous 
Sa ice bla 


are unknown, this vector would 


appear in the integrand between two shape function vectors, 


amo therefore closed form ,integration is not possible. 


When the term ECA, is included in the integrand, the 


result of integration is commonly known as the element 


Seiiliness matrix. 


In the present formulation [{k*] is called 


Ou 








the "modified e€lerent stifinese ratri. wee = KG are 
measures of element stiffness. 
Since strain energy is a scalar, the total strain 


energy in the beam is the sum of the individual element 


strain energies. Hence we have: 
n n 
Ue ee = A. (v}; cca {vi 
i=l alt 


where n is the number of elements. 
Pee OLrenvoial Bmerey.ot External Horces 
The potential energy of external forces for the eielo 


element may be written as: 


ieee ee ye 
at O at sl 
oc Heel 


The vector Sle is ~ecailed tChewelement consis ten vse acer 
Mmeroution vector and is described in detail in Appendix B=-3. 

Since the potential energy is also a scalar, the 
meral is the sum of the contributions of the individual 


pelemencs. Therefore: 


> 
SU ar 


= 
ll 
8, 
Ti ila pm) 


Ss wne Avemented Potential Enerevaerumetional 


The constant volume constraint in finite element 


i 
Ero y ae we form the augmented potential energy functional, 


n 
form Ls a ALA = a In order to maximize the load inten- 


ll! Mies 
58, 
A 
we 
V 
Feaeaeen 9 
< 
—~ 
| 
> 

Ww p 
> 
co 


n m 
a . 5 A, ae [k*] {vi ; 
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BE. THE FINITE, ELEVEN VARL ATION Le eer. 

The extremization of the potential enerey funcpwenal 
must be done globally, i.e., over the entire beam. A trans- 
formation is therefore made from the local coordinates {v}, 
meme global coordinate system {v.}. The transtormar uci. 
made by numbering the components of tv. } sequentially from 
me Origin, where the first component is the displacement 


at the origin, the second is the slope at the origin, etc. 


aos : 
o. = vi 
ale ‘IL 
= al 
Ua 
= vi = 7 
3 2 1 
ee ee ge re 
Wy ly es 


Therefore the global representation of the local "dis- 


placement" vector el becomes: 


ee ae = 1%) 


at 


oer a 
MOE 


In global coordinates the potential energy functional 


becomes: 
XC aes i ™ : 
a tv, h Uk*)| (v,) oes eb iv ae 
{=i > : i ate 
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The algebraic equaticns déseribine vanes ruler 
optimum beam shape and maximum permissible load intensity 
arise from taking the partial derivatives of the global 
augmented potential energy functional with respect to the 
state variables of individual slopes and displacements and 
also with respect to the control variables of the element 
areas and setting these equal to zero. The constraint 
equation is also imposed. 

The following set of equations, as developed in 
Appendix B-4, results. The equations must be numbered as 
meow to facilitate enforcement of boundary conditions. 


This is described later in the development. 





\ D. 
% Vv -- — = 
A a Ka oe By mee of (1) 
. Do 
a eo Oke 8 — 
ly \ 
i aces — 
i-2 aa K35 “(gei-4)* Gi/ep A 1G) Gate) 
PD,, 
ms omc |) = WG} (a eale) 
EC 
4 _ \ - 
Me Sah) GR OR 5 iE aS 
som j=l jeu 
BaD 
Xs “ee mn 1 Doin sen iy 
\ PD 
¥ om O (2n4+1) ie i 
An a 33 (j+2n-2) EC : dia 


ZT 








ho PD 
(a i a eee (ene 2a (on+2) 
ae Oo eae) EC ee 
ji 

ul 4 DXA 

2 alee aa 0 (2n+2+m) 
eal cling V(342m- 2) Ves 42m- DY ae 

Wo =) bee 

n sd 

r A, - —=0 (3n+3) 
j=. * A 

where: 

D, = a Component “Oi the flebal consistenrs eda 


is Gia Us OMn Vee vores 1¥(2nt2) 


X = strength constant 8° /6E as developed in 
Seton I. J 


Vv. = pes component of the global displacement 
vie C Cent. 


This is a set of 3nt+3 nonlinear algebraic equations 
which, when solved simultaneously, results nlectplae 
i) values for ntl nodal displacements and n+l nodal 
slopes at maximum load intensity. 
ii) the n optimum element areas 
iit) sehe Maximum load invensaLy , a 
LemoouUmGdarymecondiCiOns areulneluGed in une ereeceaine 
eapebal formulation. Since the derivatives may only be taken 
Peery respect Lo variables, a fixed displacement or slope 
boundary condition requires the removal of the equation 
associated with the variation of T* with respect to that 
moundary Condition. Ets value must be SUDSTITULeEa 1m the 
remaining equations. Boundary conditions must be consistent 
with the finite element model, i.e. they must occur at 


Megat pelnts. 
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e ) ace ‘ e ca es e 
For example, if the j state variable, ae is fixed, 


the j°" equation is eliminated from the set, and the value 
for Vv; substituted in the remedmine eauauten 

In thas way, boundary conditions are anciuded itimene 
formulation for specific prebiems. 

The set of equations with boundary conditions considered 
is the finite element model for the approximate structural 


Sorimization of a particular rectangular cross-section beam 


of uniform height and fixed volume. 
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IV. APPLICATION OF THRE SINITE Cereb 


In order to apply the model developed in the preceding 
section we must find some method for solving the system of 
nonlinear simultaneous algebraic equations. The size of the 
system for even a small number of elements requires the use 
wa Computer algorithm for solution. The only such 
algorithm available to the author was a Fortran IV-G sub- 
Routine called SUBROUTINE NLNSYS. This subroutine is 


moernded in the computer program section of This report. 


A. SUBROUTINE NLNSYS 

WamscmSUubrOuUtIMe, in theory, solves simulvarcous ly ey 
an iterative scheme developed by Brown [8], [9], [10], a 
Syscem of nonlinear algebraic equations of any order. As 
meets Cimensioned it is resGricted to, at most, 100 
equations. 
NENSYS réquires as input a vector of starting values 
Momead! | unknowns. This starting guess vector must be in 
the region necessary for convergence or the process will 
diverge. This region cannot be defined and different start- 
ing guesses may be necessary to obtain convergence. It 
should be noted that large systems of equations require 
hearge starting @uess vectors. Since the probability of 
mC OMDOMeI bE. Deime outside the region Of ,eonversence i nereadses 
with the number of components, the probability of a starting 


guess vector being within this range decreases with 
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increasing system size. (his was found (Quoc segdae sien 
that limited the size of the example problems to three 
elements. 

NLNSYS also requires an external subroutine that speci- 
fied the system of equations. RBecause the svstem of 
equations developed in Section III changes with the number 
of elements and boundary conditions of a problem, a general 
subroutine is desired that will generate the proper set of 
equations for a given problem. SUBROUTINE FCNLST, provided 


ma tCne program section, accomplishes this task. 


eee SUBROUTINE FCNLST 

SUBROUTINE FCNLST generate a set of equations which are 
evaluated at a point provided by NLNSYS. NLNSYS selects 
Only one equation at a time, thus the set of equations must 
Memoraecred in some manner. 

SUBROUTINE FCNLST sequences the equations for a given 
number of elements as they are ordered in Section III. It 
evaluates these equations at the point specified by NLNSYS. 

The equations derived from parameters that are fixed 
by boundary conditions are removed from the general set and 
NLNSYS operates on one of the equations in this reduced set. 

Problem parameters are read into the program in SUBROUTINE 


ieiisol the first time that it is called by NLNSYS. 


©. SUBROUTINE OPBEAM 
SUBROUTINE OPBEAM, as provided in the program section, 


reads the initial guess vector into the program and provides 
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the output from the program. It permits an erticiepeen ae 
of trying. multiple startine guess vectors for 4a 2ien 
problem. It starts the solution to the™ problem by sco 


NLNSYS. 


te, INSTRUCTIONS FOR PROGRAM Use 
It is first necessary to define the problem to be 
solved. We specify: 
I inks Veinene cent. wile leven (alae = SiG 
2 tne 6 ume (Gm) Ss VOU, 
3) the height (in) - HI 
4) the modulus of elasticity (lb/in*) - E 
5) the yield stress (1b/in°) - SIGYP 
Selecting the number of elements, NI, we number the dis- 
placements and slopes at element boundaries. This must be 


@eme as indicated in Figure 4. 






4 
| v,#displ Va=displ. v-=displ. a> ae 


V5=Slope v,=slope vg=slope Ve=sSlope Vo.nt+o 


Figure 4. The Finite Element Numbering Scheme. 
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We consider the boundary conditions of the problem and 
morm two vectors. One, called TPC, is 4 vector of ube: mae 
ef the components of the displacement/slape vector, V, which 
are fixed. These are ordered in increasing value. The 
second, called BC, is a vector of the values of these fixed 
eomponents. The number of boundary conditions is called 
Wee. For example, if the displacements at the left and 
ment ends of the beam are fixed at zero andthe slopes at 
these points are variable (simply supported), the vectors 


ae . 


a BC (1) 


PS ls) 0 


Be. 2) CZ SNe I Bec.) 0 


ine total number of parameters is equalage (3 xa) + 3 
and ig denoted NT. Subtracting the number of boundary con- 
Gitions, NBC, we obtain the total number of unknowns, 
me NT “NBC. This is also the size of the reduced system. 

Prom the given loading condition on the beam we form the 
Paneioat Consistent load distribution vector <D> for the number 
of elements chosen. This is denoted CLOAD in the program. 

The initial guess vector, X, is formed by guessing the 
magnitude of the unknown displacements and slopes, the 
areas of the elements, and the maximum load intensity. 

Maese must be ordered as follows: 
1) Displacements and slopes as ordered in the 
vector {v} with boundary conditions removed. 
2) Element areas starting at the origin. 


Se loxamnume kead intensity . 


oe 








We select values for MAXiT. NUM 2G, aqggee eee 
their descriptions in “SUPROUST : | ieee 
A calling program is written which dimensions X 
(the starting guess vector). This dimension is NU. The 
program calls SUBROUTINE OPREAM with numerical values for 
igo arguments except for x. A sample calling program is: 
DIMENSION x (NU) 
CALL OPBEAM(NU ,MAXIT , NUMSIG , IPRINT,X) 
SOP 
END 


*(numerical values) 


The user must provide the two statements for SUBROUTINE 
HOWLoT as described in the program. The first is the 
Same msi1on statement: 
PiMeNSlON TBC (NRG)4 BC (NBG SCLOAD( 2x12) {Gh 
Pia VCexNid? ) . AC NEE) 
Miemermensaoens are numerical values. 
ine seecona 15 the data stavement.: 
PAA NU/#/7 NBC/*/ SNL *7 SND et 7 
where * indicates the numerical values for the parameters. 
Data cards are prepared and inserted in the following 
meeecr and format : 
1) The guess vector X - 6E£12.5 
ome tone wonvsical Constants: 
PUPNC eee. Vole olGY Pee ore ee 


s)> The boundary conditions: 


a 








TBC.) ee Sate os 


TBGC2 ie. BCG2) seinen 


TBC (NBC ):.BC (NEC) eee 0. Bae 
4) The consistent oad distribupvon vec per 


CLO ADM wc alee ya5 


ine tial outous may be one of the folltowine - 
1) A singularity was generated at the output vector. 
2) The iteration limit was reached. 


Soe SOluriOn wacereacied: 


In case 1 a new starting guess should be made and the 
problem re-run. In case 2, if the number of iterations was 
preeeit< 30), try a larger number of iterations. If Gi was 
not small, try a new starting guess. Many guesses may be 
femeerrea Defore a solution is reached. 

An example problem set-up is shown in Appendix D along 


teem the results obtained for different problems. 


oe. 








V. RESULTS, CONCLUSIONS, AND RECOMMENDATIONS 





Pee hE SULTS 

Figures 22, 26, and 28 compare the three element finite 
element optimum shape with the exact optimum shape for 2 
meecitic beam under a specific loading condition.  Vable wy 
Cempares the maximum load intensities for the finite element 
designs of Figures 22-28 with those for the beam of equal 


meliume, ilength, and height but constant cross-section. 


CUB LE al 
COMPARISON OF MAXIMUM LOAD RESULTS 





Loading 











‘mo ily Supported 
Emitorm Load 







5 Oi /~in \ feeb 7 in see eel eee 





Concentrated 
imoaads at 
















x = 10 in HO Ou Kae - 194 KIP 194 
foe 2-0 in pre ae gl lela - a ONS Sallie: 191 
oa 30 in On he _ 162 Ue ike i 
mee 40. in alk KP ~ 63 KIP 183 
iene 50 in Oo Ean - 58 KIP ale 
x = 60 in B06 kale ~ 57 KIP 186 
Merancular Load 996 lb/in _ 1839 KIP 185 





Cantilever Beam 
Uniform Load 






Loi iba ee bata | Soe arene 





Maximum Load for Constant Area Beam 
Maximum Load for Exact Optimum Beam 
Maximum Load for Three Element Beam 
Poeimeresse Of Load from Case A towCase C 
pe heviation of Case C from @Gase B 


ees 
B= 
c= 
p= 
a = 
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=! +? “ys : S Tris 5 | ~ es = +} << 
It is clearly evident in Figures 22, (203 4:33 ee 


three-element design approximates the exactuep. uy ee 
Mae best criterion for the accuracy of Ghewano] 2 
mewever, iS not the beam's appearance. <Giiice OU= Foal ices 
Memimize the load intensity, the accuracy of th= <2oproxima— 


men is best measured by the closeness of the lc=2a intensity 
allowed by the finite element design to that all:iwed by the 
eeeict Optimum design. 

In those examples where exact solutions were x“xnown, even 
the coarse three-element approximation allows a >-2ximum load 
ieensity that is close to that of the exact op ti~Ume@shape. 
Table I shows that the three-element design giv=: maximum 
Weads within 5.6% of the exact design maximum fc> the simply 
Supported and cantilever beams under uniform loz: and within 
6.6% for the simply supported beam under a conc=xtrated 
ijaa at its center. 

Attempts were made to increase the number cz: elements 
meu he Simply supported beam under a uniform lces. This 
was desired in order to show that the shape wou_s converge 
to the exact solution with an increasing number of elements. 
However, difficulty was encountered in finding =n initial 


guess vector within the range of convergence. 4e to the 


{) 


limited time available and the need to check th= model's 
Zerby Lor ouner Loading and boundary Gonditt24s, the 


three-element model was used throughout. 


Sif 








Ee.) CONCLUSIONS 


The above results show that the finite ele>=nt model 
@eveloped in Section Ill, for the structural cae aie em 
of a homogeneous beam of uniform fixed height, =:d fixed 
medume, is valid. Furthermore, the results sre. thatewa 


relatively small number of elements may be usez to design 
a beam that will closely approximate the strerzth of an 
exact maximum strength design. 

Moma Orackvical matter, the consvuructon eo: —S beam ir 
the shape of the finite element solution is muth easier 
mar constructing the exact optimum shape. In the exact 
optimum design, points where the cross-sectiorz! area is 


+> 


Zero must have some material added to account Yor the shear 


effects that were neglected in the formulation. Much less, 
if any, material must be added to the finite =_sment design 
at the same points. The finite element design would have 
eitess concentrations at the element boundaric= which are 


feeepresent in the exact optimum design. 

The finite element model provides a means sf approximat- 
meee solutions to optimization problems which are impossible 
WMemcolve in a closed form by classical technizues and 
extremely difficult to solve numerically. This is due 
mainly to the algebraic nature of the finite =_ement method 
ememche diiferential equation nature of the c_assical 
Mevnods. 

The model developed was, of necessity, for the simplest 


wype Of beam and loading conditions. The apcrcach used, 
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mowever, should be valid, with some moO@hisieac emcee coer 
more complex cross-sections. It has been shown that the 
finite element techniques can be successfully applied to 


the problems of beam optimization. 


C. RECOMMENDATIONS 

Many more problems for which the developed model applies 
should be solved with various combinations of loading and 
boundary conditions for a parametric study of the beam 
optimization problem. Problems with more elements should be 
solved to show convergence to known exact solutions. 

Although SUBROUTINE NLNSYS provided solutions, the 
algorithm is very sensitive to the starting guess in rela- 
Mon to the region of convergence. More than one attempt is 
usually iEe@guimiree vo ObLaAlLnN a Peet on TO4a Biven prov lem 
Puma leorithm that does not exhibit this sensitivity is 
required to enable more efficient problem solving using 
finite element techniques. At the present time very few 
mleorithms for solving systems of nonlinear equations exist. 

Another technique that is possible would be to use some ~- 
type of direct search or steepest descent algorithm to find 
eae optimum solution to the finite element potential energy 
MmemeviOonal itself. With this type of formulation it may 
be possible to allow variable element areas. A solution 
might be more readily attained as compared with a solution 
to the system of nonlinear equations. Some work was started 


in this direction late in the investigation. The initial 
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results looked promising Dut Cimesrest iced. meee 
Bnet this approach be terminatedsin favor ot cee. ole ae 
the presented formulation. 

Models for other cross section types could be developed 
msing a similar approach to that presented in Section I1il. 
Also models with other than equal size elements or for 
Deams of more than one material would be of interest. 

In the laboratory, an experiment could be run comparing 
mae behavior of an exact optimum beam with that of its 
metre element approximations to see if the results compare 
mimes chose predicted by the models. 

The entire area of the application of finite element 
mevnods to optimization has hardly been touched. The 
finite element method is a powerful tool and more research 


mould be done in its application in this field. 
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APPENDIX A 
EXAMPLES OF BEAM STRUCTURAL OPTIMDZAL ION Fron ingy i 
USING THE VARIATIONAL FORMULATION 
i. Consider a simply supported beam of rectangular cross— 
section with constant height, h, and variable width, b(x). 
The beam is under a uniform load distribution and has 
Piven volume and length, ve and L respectively. Find the 
optimum shape, b(x) such that the load intensity 6 is a 
maximum. The beam has a modulus of elasticity E and yiels 
Serengeth Sy 
Deuce en 


mae moment of inertia of the cross-section is: 


fnus C = h@/12 and m= 1. Usime equatzon) (ll—o ) Suvi eo 


@eetcain: 


o 
Vi (bx-x ) 
7 2 
{ (Lea = ex 
@) 


a — 


Invoking equation (II-10) to obtain (II-11) the value of °% 


becomes: 


~ 
ll 
Ske 


mupstitubing the valuesmier UC, A, and 1 anvosequartion 


(TI-9) yields: 


Wy 








} ims.| 
a b 


O 


iE - 
Ape Cleese) Tobe 
O 


ihe integral can be evaluated in closed form. 


L 
if (bx-x*) = 1-76 
O 
as : 
6V. D 
MCX) ont = a Clipe ) 
which gives: 
6V p 
BCX) oot = 3 (Lx-x* ) 


The beam shape is graphically -depicted in Figure 5. The 


maximum load intensity is: 


ens. V 
ee 


max be 


This shows that, for the given problem, maximum load 
ierensity varies linearly with the height, yield stress, 


ee volume, and inversely with the cube of length. 


a Consider a simply supported beam of rectangular cross- 
section with a constant width, b, and variable height, h(x). 
The beam is under a uniform load intensity Po and has a 
given volume, Vee and length, Eee Pindethesepr 1 mumesiapes 
mex) such that the load anvensity , Poe is a maximum. The 
beam has a modulus of elasticity, E, and yield strength, 5 


y" 
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Sed bucaten ar: 


wae. moment of inertia of the eross—secrawcmen 


5 
worl ia ne 
12b 


i 





Thus C = a and n = 3. Using equation (II-8) with n = 3: 
zie 


V (Lx-x“)? 
O 


Da ey 
f (Lx=x= ) ax 
0 


La = 


Equation (II-13) gives 
O 
a= See 
= 
With A, C, and n above, equation (II-9) becomes: 
Sue i 


PO 7 3b ) L re 2 
E (Lx-x ) ax 
@ 


Mmreewintesral may be evaluated and is: 


L ie 2 
if Gime, ass a uae 
0 8 
wes : 
uy e 
oe te 
BV | (LX-x ) 2.55 V_, (Lx-x ) 
Va. = Sa 
ON nL Le 


Selvineg for h(x): 


255 Vo (Lx-x*) | 


bLe 


Jalna) 
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h(x) is grathically depicted in Figure 6. The maximum per- 


fess ible load intensity ise 





64s ve 2.178 vo 
MO ae 7 5 r Zi hy : 
Seo bL 


Here Oe varies linearly with yield stress, para- 
bolically with volume and inversely with width and the 


fourth power of length. 


fee consider a simply supported beam of circular cross=— 
peeuwlon under a uniform load distribution. It has a given 
momtume and length, ue and L respectively. Find the optimum 
Pedius, r(x) such that the load intensity is a maximum. 
mee beam has a modulus of elasticity, BE, and yield strengin, 
oa 
Y 

Soluu Lon : 


ijme moment of inertia of a circular cross-section is: 


n 


= 


u 2 

© Picton wilt 2 = 

I = Ses one (ir <) aes 
Thus C = i= and n= 2. With n = 2, equation (II-8) becomes 


D 273 
Vo (Lx-x ) 
ie , 2/3 

i Ob ere): dx 

0 


(oe = 


Mewacion (lI-12) gives the value of XA’ for a circular ecross— 


section as: 


yu 


With A, C, and n as above, eauation (2-9) becomes: 


é O/ 3 
Bo EVOL | ee 
O — L O/ 3 By) 
ev E Gi ax| 
0 


The integral may be evaluated in closed form. However it 
involves the Gamma Function which must be approximated. 


Mae integral is, then: 


L Bes 
ai Olea) ax = 20 7o (1) (73 
@ 


wus : 


Cis 
Vf Lx-x ) 


7 


Cig) (ME 


Solving A = mr? 


POX) oot 





The beam shape is depicted in Figure 7. The maximum 


ieee intensity is given by: 





PO : 
oS oa 
ete) a8 linéar 2) Victor eee ewer ooo omeee: Cuma 


max 


7/2 power of volume and inversely proportional to the 7/2 
power of length. 
Reference [ll] gives optimum beam shapes for example 


problems 1 and 2 derived 1n a different manner. Tne shapes 
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unitorm load cistrimition 


Pe TS intensity P. 


i 





B O a i 
eh iG) =a ee pe at 
bI, 





Figure 6. Optimum Peam Shave for a Pectangular Cross- 
sect von Peano! ni form’ iar wmeer 
Uniterr Load. 
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ee 


Circular cross c<eerten 
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eae 


1/7 ae 
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LPP awn 
CALCULATIONS FOR TPH PINITE aaa ee ae 
FOR BEAM STRUCEURAL OPT IMI ZAaT ro} 
ime ohe Cubic Disolacement Funetion 
We wish to find @ cubic displacemenu  funcusom vecten 


Sex, )> Swen (tia 


v(x, ) = <N(x, )? 


Goa Ga 
mM Mm HK 4 


where v(x, ) is a cubic polynomial. 
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Figure 8. The Deflected Element. 
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We let v(x, ) = ax; + Bxé + YX, + 6 over the interval 


we X, 


; <A. Boundary Seon dere iese re: 


Bie) 








(1) v(0) = v 


(2) v'(0) = 


ih 
(3) ov Gl V5 
ee 0. 
miipecituting boundary condition (1) in v(x, ) gives 6 =v 


1: 
Taking the first derivative of v(x, ) with respect to 


Xs yields: 


Ps 
t = 
V (xX, ) 30x, + 28x. + ¥ 


Boundary condition (2) in v'(x, ) yields: 


Seestituting the values for y and 6 into v, and v’ and 


mieiines boundary cenditions 43) and (4) yields: 


a 
A A A A 
: 3V5 e0. 3V 5 8 ., 
amiees (0 Sie 
A - A A 


Q 
i 
/\ 
I 
= 
I 
ne 
> = 
NO 
V/ 
= 
< 
=, 
—_ 


= 
H 


<0. 1 Oe {vl}, 


Or 
I 


Gi’ 0 0 sane. {v}. 


imhcwin VeCLOr MOCavLIon: 


ipl 


























wy < 2 = a) = | 
Ox 3x- toe 2. ee = ea 
ae < 5 7 Dt a ee 
A A ‘A - ; A? Am 
dx) 
meee |e 1 
Ac A al 
Therefore the shape function vector <N(X, )> Lee 
ae 2h x2 24° 2x2 Ko 
<N(x,)> = Se 
A A A 
ba. 3 5 
x 
o 


Ke 
wl » 
A A 
ieee second derivative of this vector with respect to Xs nish 


Bee /l te Se ae 
all Ne ae Ac We ne Ne A 


2. The Modified Element Stiffness Matrix 


Bh aes 





> 


Mae modified element stifiness matrix 2s defined 2s: 
- ie 

fee |) = f SNe esc yl dx, 

Tera 0 4x1 1x4 


Taking the indicated vector product and integrating yields: 


[#] = 12/A° 6/A° 1 2/A? ie 
6/A° N/) - 6/A° o/h 
-12/A>_-6 7A" 12/A> 6 A 
6/A° BSA SAE OYA 
fee Phe 6 Consistent Load Distributions vec ver 


Assume that the load per unit length over the i 2 element 


may be written as an oe Orde repel Motmras 
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mmere Oo, B, Y;«..,nN are constants Whiten sreugeve video eeege 


leqewn Values of load intensity ab rtl loacvions a2 





O< x, <A. These known values are represeitve¢ asm ocr mcne 
of the maximum load intensity, ee over the beam (0<x<L) 
a ecp, pee 
olen as a = a 5 where 0 < 2 aay alc 
O = 
P(x.) = P_[a + 8x. + vxe + + nx? ] 
iL O ae ae aa 


The load on an increment ax, of the bh element is: 


= a r 

dP(x, ) = P jlo uO eee nx, ] dx, 
The work done on the ae element by this incremental force 
mS 

dw, = dP (x, ) v(x.) 
uc. : 

v(x, ) - <N(X, )? Hae 

pe! Lm 

mous : 


_ 2 ie 
dw. = P fatBx tyx,t...+nx, | <N(x, )> {vl}. qx, 


The work done on the qe element is: 


A 
Ne if dw. 
0 
= fa Po leases xo + + aoe <N(x eo Gx 
oa nes ee he i 4 Aa 


a 








7 


betmanc Wendy Pam ea ani ps | be denoted by oy ana 
meting that i and {v}. are independent of x t- = above 


expression becomes: 


A 


as R ; eC) <N(X, )> dx, {v}, 


The integral is defined as the element consic=ent load 


a 
= 


Ses tpribution vector ay 


A 
<D>, = : (x, ) <N(x, )> dx, 
1x4 pe 
The work done on the oe element may then be wriz7en: 
Cla ee eek Oe 


The work done on the beam is 


n n 
W= YY w. =P i D> ae 
: il One ut al 
fee i=1 
It may be noted there is no requirement for the =_ement 


Mmeead distribution polynomial, Ae oa GO be the S=se for "each 
element. On the contrary, each element may have = different 
ieea distribution polynomial, P(x, ), aiid, sce dikes 1 Ue ater 
may even be discontinuous at element boundaries. The only 
requirement is that the ratio of P(local)/P, be =xmown at 
met DOIntS over the element in order to define Tze eae order 
Seesouribution polynomial. 

We take the following case as an example: 


Consider the first three elements of an - element 


peam. Given that the distrapueiens of load on —-=-;cultc =. 
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¢e, and 3 are cubic, linear wane see 2s - respectively, and 
also given the values of eee where 2 = £9085) (2 
we seek an expression for the work done om the fires circ] 


elements of the beam. Figure 9 shows 4 graphical represen- 


meacion of the problem. 





Figure, 9. Multiple Loadans Uy yeecuon aahecin 


The values id, / Ee Chnrouecin CD ee define the cubic distri- 
meeion function 92(x,) over the first element. The values 
PLP. and PAP. define the linear distribution function 
Q, (x5) over the second element. The values of PL/P 
ao OuUrh El define the quadratic distribution function 
5 (X) over the third element. These functions in turn de- 
mene consistent load distribution vectors <D.>, <D.>, and 


a 2 


<D.> for the elements. 


e 
A 
<D,> = : 920%, ) <N(x,)> dx 
A 
<D,> = ; Q, (x5) <N(x, )? dx. 
A 
<D.> = : Qo (x3) <N (x3 )> dx. 


5D 








m4 bh =| = e P <p eg ; Je = . 
fhe work done on the first three ele on. 


is (Pree ee Plo s Ue + <D>, {vet +02 > el 


2 Ae 


4, The Concentrated Load 

if the load 1s concentraved Vice d1s7 10 ooe a 
wonsistent load distribution vector is formed Miwa emer 
Meeailel to that of the disuributed load. 

We apply a concentrated VGad 67 sinvemsery BECO <p) 
at a distance aA(0 < a < 1) from the origin of the th 


element. 


pk 





Pigure 10. A General Element with Concenrravecsload: 


mae work done on this element is: 


Ws = BP vad) 


EU: 


wah). = <N on fv}, 
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Mherefore: 


Cheat ale <N(aA)> {vl}, 


If we define the element consistent load distribution 


vector as: 
<D>, mae cei C ow iy) = 


th 


we obtain for the work done on the i element by the load 


ees 
O 
= > 
We - <D>, tv} 
Note that the expression is the same as that for the dis- 
meapuved load. However, the units for Pe and “Dae are 
different - ae ToL: bomwOr LLOrCenViCesPorce speci unin: 
length. 
imi muller ple “Concentrated Woads are applied to lone 
element, the element consistent load vector is the sum of 


Chose obtained by considering each load separately. In 


symbolic form: 


Wine individual load vectors are: 


<D)>, = 6, <N(a,4)> 
<D5>; = 6. <N(a,A)> 
SS = ee es 

ia B (0. ) 


the ebement consistent load vector is: 


<> ST ESIN Ye) aga 
1 = EB; <N(a5A) 


ai 








mM Toate Of = PE ay 2 7a ‘Cc : 7s 
om The Finite EPblereants Sriatlitone | -ope 6 .- 
se ge gg ee ee em eee ee 


9 


We have as the g€lobal augrented potential energy 


mune t Oneal: 
n ial 
Tf ee Ge [k*] {v.} - P my <D>, {v.} 
g oy sl a 1 a ores 1 1 
j=] j=] 
i 
= AA. 
i=l 


where n is the number of elements. 
Define the following as the global consistent load 


meee r 1 OULLON Vector <D> 


ia 2nte) 
<p> = <pl. pe. p3 + pt. pt + pe pt. + p2. p3. nis 
7 Pa Pa Be Ps poco ekg eg 
= <pD_, D., D i. 


q? Pos Pas+++>Dan47> Ponto” 


<D> is formed by assembling the element distribution vectors 


as indicated. Thus the global work term becomes: 


le 
O 


[ape 


SG - D> v7 
<D, tv. } Be <D et 


se 1x(2nt2) (2nt2)xl 


aL 


The augmented global potential energy functional that will 


be operated on is: 


iQ) 
a = : A. {v,} [k*] {v,]} Be i ed va AA 


lec SS 
‘T> 


Taking the partial derivative of T* with respect to each 


component V, ier) sandesetting rescues |. cee rome coe 


oT # y 
= eye Chee: ce 








dV 5 1 = ia) puke O 2 
aT# \ 
= = ECA 4. k& ae 
- + 
8V 44) (Ase) eo 353. ° (jt+i-4) 
a = 
x 
PEA 72) 42, “Tg V+4-2) Giza 
Ho We en 
ats : 
——— = ECA kf 
OV. (ize 2) fe a hi; (jt+i-4) 
4 a 
% 
eA a. , “Aj thay > ae 
i = Gyo. on 
oe 4 
aa can ee erage 
°V (antl) Mo jey 3d © G#2n-2)- Py Deonagy = 
oT y 
ame be Kee Wore mi = 


Mis £ives 2nte equations. 
We emexy Take Pine partial derivatives, of TS with respect 


pee one A,'S and set them equal to Zero: 


oT ul ut 


A 
: 2 shee : aa 
OA. jock a 1j (kt2i=2) WC jes Oe ane 


yee ae i 


This gives n equations. 
At this stage we have 2nt2 equations in 3nt+3 unknowns 
Mence anvaadivional equation is necessary. Ths =equavicon 


ino ende Volume COnStLraint : 


De 











Therefore we have a total of 3nt3 equations in the 32n+3 


unknowns: 


2n+2 components of the vector {v} 


igh Buc Ses Ie 
% 
1M ike cay eit Iherskel) setegver sya) cay Ee 


We therefore have a mathematically tractable problem. 
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ja ecules oo BJAG, 


BAAMPLES©OF CONSDSRE WP LOsAD Dist ee UO VCR 


1. As our first example we consider a three element problem 


maaer 2 uniform to0ad> distribution: 





Figure 11, A Three Element Beam, Uniform Load. 


Olserving the general element of Figure 12, the work done 


Samo nis element is: 


A 
We : P <N(x, )> dx, {tv}. 





Figure le. <A General Element, Uniform Load. 
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which may be written: 


A 
= Pp i . 
W. P 6 / SCO: dx, oar 


From our previous definition in Appendix B, the element 


@oncistent load distribution vector is: 


A 
<D> - ; <N(x, )? dx, 


Evaluating the integral where INO i Sclsordle ta me cmt 


Appendix B: 
A 2 2 
ee ee 0 es 
Z MCS es Sa res oe Fay > 
aus 
2 2 
i ee 8 
ar oreo sae (o= 1) 


Peeevying this to the three element problem and noting 
i@gaw the load intensity is the same for each element, we 


obtain: 


<D?, = <Des = peas =<3 


We assemble these element vectors to obtain the global con- 


Sistent load distribution vector <D> 


<D> =< : > Cee?) 


The work done on the beam is: 
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<< | 


eee 

ite. 

2 2 V 

A A A A MN 
eae TD? NO oui ee 2 me Ve 
"6 

af 

Vg 


Oo. AS a second example let us take the following problem: 





Vig 


Figure 13. Three Element Beam, Multiple Uniform Loads. 


Mmpe element consistent load distribution vector for a uni= 


memm load from equation (C-1) is: 


ae e 


*) 


rol 
Ne 
ed a 
NO 


AN 

S, 

V 

It 

N 
nol > 
| > 
NO 


In this problem, however, the load intensity is not the 
Same on each element. Therefore the element consistent load 


Meso UrIOULION Vectors are: 
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tH] 

) 
Pole 
ed 


O a 
One: 
2 2 
A A A A 
2 2 
Bees So ei aS 
SO TD DS Para 


Assembling these vectors we obtain for the global consistent 
iad distribution vector; 


A Ac 3h Ac A one 


an A2 3A Ae A 
oT BI 9) Te Wa Be he aie 


bee Let us now lock at the following problem with a concen— 


trated load: 





Figure 14. Three Element Beam, One Concentrated Load. 
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Figure 15. General Element, Concentrated Load. 


MmeremcONStsred) HOnd Geet ri burl On VeCcuorerer biad sircameiie ayes. 
<D>, =a Can 

But : 
<N(aA)> = <(20°-30°+1), (a2A-20°atad), (-20°+3a°), 

GENS > 

aus : 

<D>, = <(20°-30°41), ican daceasen) ae iar aeons ne 
[(a°A)(a-1)) > (O-4) 


In the specified problem there are no loads on elements 2 


orm 3 thus: 


Des = <D.> = sic 0. OG. 0e 


The e lemme consistent load vector for element 1 is there- 


fore: 


NO 
—~] 


i 3A = S, 3A 
<D>, = i = Ble 6H? 


ao ai > (C-Ha) 


8) 
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1? ar and <D.> we obtain the ¢ilobal con 


sistent load distribution veetor <D-. 


Assembling <D> 


<D> = < 25, Jf, &, - A, 0, 0, 0, 0> (C-5) 


i) Let us take exemple 2 andeaddwan edeit 1onel sloadeon 


element 1 of intensity A a distance : from the left end: 





(a 


yoV 


Figure 16. Three Element Beam, Two Concentrated 
Loads on Element One. 


Tne element consistent load distribution vector for a 


concentrated load ee a Gietencenk, lt rome themerd oats 


ear _ af oho ZA 
“Daa = AGS 6? 2 os 64 ? (C-6) 


In this example <D>. and <D>. are the same as in example 3. 


@ 3 


iat is ee = See = ~“UnN0nOLo>. Ihe Consistent Jo0edsdis— 


tribution vector for element 1 is the sum: 


<D>, = <D(f)> + <D(q=)> 
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=) 
< 
C) 


hy 2f 


a 


2: 
5 D2 eels il ’ le = Papas 


Assembling the element vectors we obtain for the global. 


vector: 


A 
<p> =) <1, 2, ie =, GER Ces On 10.22 (C=7) 


5. For this example let us take the following beam with 


eoncentrated loads as shown: 





Figure 17. Three Element Beam, Concentrated 
Loads on Elements one and Three. 


From equation (C-6) we have for <D>,: 


eae ee 3A 
acig” mg) GHEE rine 


Element 2 has no loads thus: 


SP acs On On Oe 


mawacion (C-4a) gives the element consistent load vector 


for a loading such as is on element 3 but of intensity ae 
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The element consistent load distribution vector for element 


3 is therefore one—-hadf of vhav vee jcece= i eee 


mee, OA 
<D>3 = < gms ipo By? ~ 438 > 


Assembling these to obtain the global consistent load 


suri bULLonm Vector: 


5 34 27 9A (0-8) 


ee oo” 2 
ae 0) 32: 35> - Bl» by» tow bee - Tae? 


2g 


wees a final example let us take a three element beam 


Pager a triangular load as shown: 








\ : a , | . 1 2 pe prea ; —s 
Ak A ele 
vo 


Figure 18. Three Element Beam, Triangular Load. 


We first examine a general element with a basic 


triangular load of maximum intensity P as shown: 
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Figure 19. General Element, Triangular Load. 


Hoes distribution may be written: 
oe 
ait 


aoe 


The element consistent load distribution vector by 


@efrinition is: 
pe <N (xn) a 
, a al 


Evaluating the integral we obtain for the basic element 
Sensescvcent load distribution for a triangular toad; 


P(x, ) = PX. : 


> (C-9) 





A —re— A tle — 4 ——r| 
Pigure 20. Three Element Beam, Superposition of 
| Dia lCrieena rriangcular Loads 
¥oV 
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We see that element 1. has a bdasie triangular o4dedi er. 
bution of intensity Po/3- Element 2 has a uniform load 
Secor i bution of mavens y ae ae. with the triangular distri- 
moa On Of element il Superpesed. Element 35 Similariy has 
ais SUpDeErDpOSition but, the umiform toad anvensity is 

a / 3. 


Peon C—o MancGupmeere lapicns: == eee We cot alin. 


—] 
> 


P O 2 
2S € al <p> vriang — A. 


A 
lo 3 zt 20° 90° 60- 


I 
V 


NO 
CO 


From the result for <D>, above, and from (C-1) with 


mocensity ee EhesUnit orm Portion Of BES se 








eppumiform 1, A AS A =Ae | 
2 Be 1 oe 
mus}, 
paumnitomn a, Die i ne © 
Ce 628 Canale 
And: 
<D>, 2 <p> eriang fh Desa 


coe 2e8 TE as ae 
2 Con? TO? COp A Se 


The third element consistent load distribution vector is 


developed in a similar manner: 


A Oo o 
gal it Og (lee 24 A AS A | AS 
<)> = 2 < 72 [ox D> 1S > 
Ns ee oui 
Sees as 
And: 
SD = al F: i alae 
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ine: 


Assembling <D>, ; <D>,; and <D? 2 we obtain: 


bo Ae 2k a oe 13a , (C-10) 
AO? Sie a ys? 2° U5? ae a 180 re 


(ue 








ARP EN DL Koel? 


PINLTR SER EME Ne Orval eae rie 


All of the following examples consider the approximate 
Svructural optimigation of a rectangular cross-secrt von beam 
of uniform height and fixed volume. The physical parameters 
are: 

a) beam length = 120 in 
D) beam heleht = 6.4n 
c) beam volume = 2200 in? 
gd) modulus of elasticity = 30x10° 1b/in® 
e) yield strength = 50x10° 1b/in*® 
A three element approximation is used throughout. 
im Te pleco Ua Ve spr Oo wen ysel=Up 
As a representative problem consider a simply sup- 


morcved beam under a uniform load. The beam is depicted in 


eure Cl. NI = 3. 








a 
cmnee a =, Sees nee ene ree, ea 
th. + a zs 
— eee Se nes : ; ; j / ; 
: YATT 
HQ iy —!} 0) inte 0 ina 
V 
% V2 Ve a 
2 vy VE Q 


Pigure 21. Three Element Problem, Uniform Load. 


ee 








Since the beam is siroly supported the boundary condi- 


t10ns are: 


HO 0 
V7 = Q 
Therefore: 
IBC(1) = 1 BOGE) =" 0.70 
TBC?) = 7% BG C2) == 100 
NBC = 2 


The total number of parameters are: 
NT = 3 x 3+ 32= 12 

Me coval number of unknowns are: 
NU — eee 0 


teem Appendix C, the consistent load distribution vector for 
a uniform load on three elements is: 


: an 
2 2 


A 
Seat ay ee Os 


onneyWD oa 


< 


Mj 
rol 


Eaarce A = L/3 


Homan eo lbOn Ds becomes. 
<P Olbs ceaas So tO.) Cems O20 lo ee me 


The components of the x vector are: 


= < : 
X Xb > Vas Vivo ve Vee Ve > Ay» Ay Aas Py 
We guess values for the x vector: 
= 92 A, = 15 in* P. = 760 lb/in 
= aL 75) A, =. 26 in’ 
_ ee 
Vy = On: A. = 15. an 
mie = .9 in 
We = -.01 Vp = -,03 


ies 








foe Ta vrm Sey TtS Top We ee a fe] 
LOR a ae EUPRSIG and 122. Wey ae 


require 50 iterations maxinun. 
MAXIT = 50 
me cesire four Silent iicanteateirt.. 
NUMSIG = 4 
We desire the x value printed at each iteration. 
TP Ne) 
Me miOow Can utilize the program. The calling program is: 


DIMENSION X(10) 
CAVE OPER AMIGO. SOO ) 
Sa Or 


END 


Meminsert the tollowing cards in their proper locations in 


FCNLST. 


Die Nomen aC C2 ABCC?) CLOAD( Ce G. lial. 
Coe? 
DADA Ni Oya C/ 27 ste ey-aNihy ey, 


The Aiea data Card 1S’ Umemwtirst Sixeeemponenks of the 
Paevecuor in Fl2.5 format. The second is the last four 
mareOnents of the x vector wim Hl2.5 format. The third card 
iemcie dast Of piv sical *parmamevers PLENG.E ile vObesIGYP an 
tie. > foOrmab. Wie fOurthwecard is the first “component ef 
WBC and the Corresponding scomponent of BO (il0 Bie 75). 

The fifth card is the same as the fourth but for the second 


Compcnents of the vectors. 


TH 








The last two cards Contein the e1eay Compo tenes roe at 
eonsistent load distributionmevecctor 1a 521225 ore 


data deck is thus: 


O.. 3B=1 OOO O11 R= 0.9E0 S18) AS ~0.3E-1 
ie oie! ey olde easy Jk eo 2 
ee ok Bley 6.080 Pagel 5) (Ole 

i O7..0B 0 

2 G0b0 
2 OBL eee So 5) 2), Os 4.0E1 Cr Oia 4.0E1 0.OEO 
2 0dsil eee 35) 2 


This program was run and the following solution was 


leeehurned aiter four iterations. 


Vy = .0310 ho = 5.0L a Pa = 20h cn ley ene 
Vz = .9566 in A, a 2s in’ 

oy 2 Oui A, = 15.01 in 

Ve = .9566 in 

Ve = -.0111 

Ve = ~.0310 


We are interested in mainly the areas and maximum load 
results. If the areas are divided by the height (6 in) we 
Seeain the width. Since the beam is symmetric about its 
@eamcerline we may plot one half of the beam Co depict its 
pape. Figure 22 is a plot of the three finite element 
approximation with the exact solution superimposed. 

poe OlLNerre amp te pee lems 

The problems of a simply supported beam under con- 


Gentrated loads at various locations were also solved with 


ie 








the computer program. The results for these problems are 
shown in Figures 23, 24, 25, and 2. Corresponding loads 
to the right of the beam midpoint were checked and the beam 
shapes obtained were the inverse of those Shown, i@e., the 
half-width values for the end elements were interchanged. 

The solution of a simply supported beam under a tri- 
mepular load distribution is shown im Figure 2/7; 

fanealiy the problem of @ Canui lever beam Undersea unicorn 
ieac was soived. This is shown in Figure 20. 

Optimum beam shapes for the simply supported beams under 
Mmmeominitorm load and 2) a concentraved load at the center 
and also for the cantilever beam under a uniform load are 
shown in [11]. These exact optimum shapes are superimposed 


fiomecomarison in Figures 22, 26, and 28. 
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77 as 
6 in 
/\ VOLUME = 2200 in3 & 
FE = Cee lb/in* ohd 
os = SOLO. We Hie 
| yy i ‘ X 
0.0 2n uo ~~ 60 20 LA0 120 (in) 
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ee oO Peact Op paren 


, val lye te 
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nea 
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= Design Ne 
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VEO UVEOAD , HIME LEME Ne SOU Seed 2 lo yan 


A GOT) el PACUG SOIT Tae. & a euseare «4 ee (on Loven 
MAXIMUM LOAD. UNIPORM CROSS-SECTION. .. 2. SAMO ey al 
PIigure 22. @harecseterche  pLsrum Pesien for a 


Simply supported Pear under 
Uniform Load. 
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ae ea be e-wiP rH in) 


er 


DO 
nN 


PO 
<a 


po 
Com, 


ich 
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) 
WI 


4 
ye = - 


pw X 
20 30 4O 50 &4n 20 m0 (in) 





Pmax = 194 KTP 





PeOadomaveapoints | Uhroueh 6 apply in Picurcemese= 2n- 


The beams 


Pipure 23. 


Shown in these figures have the parameters: 
Vo. = €208 in 


O 
z Ae they ea- 
See I ane 


Ss) 
v 


Three Flement Optimum Design for a 
simply Supported Peam under.a Con- 
centrated@icad Applicd arson i 
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PEAM HALF-WIDTH (in) 
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Pmax = 105 KIP 






i ee 


Pmax = 76 KIP 
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Figure 24. Three Flement Optimum Peam Designs for 
Simply Supported Peams under Concen- 
Uraeved beads .at Poincs 2 “ances. 
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Figure 25. Three Element Optimum Ream Designs ior 


Simply Supvorlted Peams umder Conecen— 
{rape VMorcs fpolaed ater ouncsu omer >. 
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Peact. Optirun 

Shane 

Prax (exact) = 
61 KIP 






/ Pmax = 5/7 KIP \ 
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i, Sauer ni pee iaaoee = ee 


mm 


Figure 26. Three Flement Optimum Ream Design for a 


SIMD ly Weeertend eam ander a Concen— 
trated Load Applied at Point A. 
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PRAvonALP=Wi DTH Cin) 
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Fares 2 7. 
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P ymax =. ee Oo Ain 


Three Flement Optimum Ream Design for a 
Simply Suoported Ream under a Triangular 
Load Disticieucione 
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0 20 4O G0. RO 100 120 (in) 
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BRxact Ootimum Ream Shape 
We Pores (exact) = 382 lb/in 
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P ymax = Some ta 
. 36 


eo omenfee 


Figure 28. Three Flement Optimum Beam Design for a 
Cantilever Peam under a Uniform Load 
Daiser ip Utica. 
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